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ABSTRACT 

We present mid-infrared spectra of 5 submillimeter galaxies at z = 0.65 — 2.38 taken with the Spitzer 
Space Telescope. Four of these sources, at z < 1.5, have strong PAH features and their composite 
spectrum is well fitted by an M82-like spectrum with an additional power-law component consistent 
with that expected from AGN activity. Based on a comparison of the 7.7-/im PAH equivalent width 
and the PAH-to- infrared luminosity ratio of these galaxies with local templates, we conclude that these 
galaxies host both star-formation and AGN activity, with star-formation dominating the bolometric 
luminosity The source at z = 2.38 displays a Mrk 231-type broad feature at restframe ^8/im that 
does not conform to the typical 7.7/8.6 /im PAH complex in starburst galaxies, suggesting a more 
substantial AGN contribution. 

Subject headings: infrared: galaxies— galaxies: starburst — galaxies: AGN — technique: spectroscopic 



1. INTRODUCTION 



Deep submillimeter-wa ve surve ys (ISmail et al.l 



1997 



Barger. Cowie fc Sanders] 119991: lEales et all 1 1999 
Cowie et all 120021: IScott et alj l2002t iBorvs et al.l 12003 
Webb et al.l l2003bl ) have uncovered a population 
of ultra- luminous infrared (I R) galaxies (U LIRGs; 
L T r>2 - 5 x 10 12 L Q ) at z~2 (jBlain et al.ll2002t ). This 
observationally-defined population of submillimeter 
galaxies (SMGs) coincides with the epoch of peak global 
star formation and quasar activity, with a significant 
contribution to the global star formati on rate density 
at z = 2 — 3 not traced in the UV (jChapman et alJ 
120051 C05). Highly obscured by their dust content, 
the astrophysics of SMGs and the nature of their 
power source remain a ch allenge to address at o ptical 
and near- IR wavelen gths (jChapman et al.l 120031 : C05; 
ISwinbank et al.ll2004l ). Deep X-ray studies suggest that 
~ 28 — 50% of SMGs host an active galactic nucleus 
(AGN), although at face value it appears that the AGN 
does not dominate the bolometric luminosity and that 
powerful starbursts (SB) contribute more signifi cantly 
to the total energy output (|Alexander et al.l[2005l . A05). 

Less hindered by obscuration than shorter wave- 
lengths, the mid-IR region boasts a number of spec- 
tral features, including: emission from Polycyclic Aro- 
matic Hydrocarbons (PAHs) (e.g. restframe 6.2, 7. 7, 8.6, 
11.3 and 12.7/im), associated with star form ation dHeloul 



11999ft and typically absent in powerful AGN (|Voitlll992fl : 
silicate absorption at 9.7 and 18 /urn, which gives a mea- 
sure of the obscuration by silicate dust grains along the 
line of sight to a small hot dust continuum source; and 
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a hot dust continuum (A < 10 /im), likely to be domi- 
nated by an AGN. The strength of these features have 
been used in mid-IR surveys with the Infrared Space 
Observatory (iGenzel et aLlll998l: [Rigopoulou et alJH999t 
iTran et ahfeOQlULaurent et al.ll2000ft to estimate the rel- 
ative contribut ions of SB and AGN for the brightest local 
galaxies (e.g. iRigopoulou et al.lll999ft . The spectra of 
the high-redshift population remained unexplored in the 
mid-IR, until the advent of the Spitzer Space Telescope 
and the unprec edented sensitivity of the Infrared Spec - 
trograp h (IRS: IHouck et all l2004f) iLutz et all (l2005f) 
CL05L lYan etall (I2005F7Y05L iHouck et alJ (I2005T). 
iDesai et alJ ()2006ft and lWeedman et al.l (|2006ft have been 
among the pioneers in using IRS to study the mid-IR 
spectra of luminous sources at z ~ 1 — 3, extending to 
higher redshifts the analysis that was previously only ac- 
cessible for nearby galaxies. 

We have an IRS program to study the range of mid- 
IR properties of a sample of 24 high-z SMGs with 
S24/jm ^ 0.4 mJy, using the radio-identified sample 
with spectroscopic redshifts, compiled by C05. Here 
we present IRS spectra of the first five targets ob- 
served (SMM J221733 +001120, SMM J163659 +405728, 
SMMJ030228 +000654, SMMJ163639 +405636, and 
SMMJ163650 +405735): four are at lower-z, with z = 
0.65 — 1.5, and one is at z — 2.38. The low-z targets 
cover wavelengths longwards of 10 /im and give insight 
into the longer mid-IR emission from SMGs; the full sam- 
ple is more focused on z ~ 2 SMGs and hence probes 
shorter restframe wavelengths. This preliminary sam- 
ple is otherwise representative of the SMG population, 
in terms of bolometric luminosity, dust temperature and 
submillimeter-to-radio flux ratio (C05). 

2. OBSERVATIONS AND REDUCTION 

We observed each target using the low resolution Long- 
Low (LL) observing mode of IRS (R ~ 57 — 126) at two 
different nod positions for 30 cycles of 120s each. We 
cover restframe emission longwards of 6 /im to probe for 
PAH emission at 6.2, 7.7, 8.6 and 11.3 /im and for silicate 
absorption centered at 9.7 /im. The data were obtained 
between December 2005 and March 2006. 
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Fig. 1. — ID Spitzer IRS spectra for 5 SMGs and the composite spectrum of the PAH sample. For the 5 individual spectra, the lower 
curve represents the unsmoothed spectrum, with the first order (LL1: A f, s = 19.5 — 38 fim) and second order (LL2: A jj S = 14 — 21.3 /im) 
of the low-resolution mode shown in dotted and dashed lines, respectively. The upper curve shows the spectrum smoothed by 3 pixels and 
offset in flux for clarity. The various wavelengths of PAH emission features are indicated. We show the smoothed version of the composite 
spectrum for the PAH sample, together with the ISO SWS spectra of M82 (dashed line), smoothed to the resolution of IRS and normalized 
to the 7.7-fira PAH feature. The excess i n the S MG composite, when fitted by M82, is consistent with an additional power-law component 
emission from an AGN (dotted line; see i|3.1.2t . 



The data were processed using the Spitzer IRS S13 
pipeline 8 , which includes saturation flagging, dark sub- 
traction, linearity correction, ramp correction and flat- 
fielding. With a slit size of ~ 10.5 x 168", IRS does not 
resolve the SMGs spatially, and the targets were treated 
as point sources throughout the data reduction and anal- 
ysis. We performed additional reduction of the 2D spec- 
tra using IRSCLEAN 9 to remove rogue pixels, and relied 
on differencing between the nod positions to subtract the 
residual background. We used the SPitzer IRS Custom 
Extraction (SPICE) 10 software to optimally extract flux- 
calibrated ID spectra, by taking a weighted average of 
profile- normalized flux at each wavelength to increase the 
S/N of these faint sources. 

3. RESULTS AND DISCUSSION 

The mid-IR spectra of SMMJ221733, SMMJ163659, 
SMM J030228 and SMMJ163639 show moderate to 
strong PAH features (see Fig. 1), and we refer to these 
targets collectively as the PAH sample. Detection of 
PAH emission is assumed to indicate the presence of SB 
activity. At most a very shallow dip is present around 
9.7 /im in the spectra, indicating little silicate absorption. 

Our highest-z source, SMMJ163650, is somewhat dif- 
ferent to the other targets, with a broad feature at 
restframe ~8/zm, unlike the typical blended PAH com- 
plex of the 7.7- and 8.6-/mi features found in SB galax- 
ies. It is more remin iscent of the spectrum of Mrk231 
(jArmus et al.l l2006bl A06b) , which features an unab- 
sorbed continuum between absorption from silicates at 
longer wavelengths and hydrocarbons at short er ones 
(|Spoon et all 12004 S04; IWeedman etaD 12006( 1 . This 
similarity suggests that SMM J163650 has more substan- 

8 http:/ /ssc. spitzer. caltech.edu/irs/dh/ 

9 http: / /ssc. spitzer. caltech.edu/archanaly/contributed/irsclean 

10 http:/ /ssc. spitzer. caltech.edu/postbcd/spice. html 



tial AGN-activity than the SMGs in the PAH sample, 
as expected from the presence of a strong CIV (A1549) 
feature at restframe UV (C05) and a broad Ha com - 
ponent (~ 1753 ± 238 km s _1 ; ISwinbank et al.l 12004ft . 
both revealing the unambiguous presence of an AGN. We 
discuss the properties of this source in more detail in a 
subsequent paper discussing the full sample (Menendez- 
Delmestre et al. in prep.) and concentrate here on the 
median properties of the SMGs with clear PAH emission. 

To get an insight into the physics inherent to SMGs 
in our PAH sample, we compare their spectra with ex- 
tensively studied local te mplates: the AGN s Mrk231 
(A06b ) and NGC 1068 (ISturm et a p [200l . the SB 
M82 (jForster Schreiber et all 120031) and the U LIRGs 
Arp220 (A06b) and NGC 6240~ (|Armus et al.l l2006al 
A06a) . Arp 22 has been a favorite template for high- 
redshift SMGs ((Pope et alJl2006tlKovacs et al.ll2006h : it 
has strong PAH features, indicative of SB activity, and 
a steep mid-IR continuum due to a heavily obscured nu- 
clear component inferred to be responsible for the bulk 
of the IR luminosity (S04). AGNs have been identified 
in both merging components o f NGC 6240 but SB d om- 
inates the total IR luminosity (jKomossa etalll200ll . 

A qualitative comparison of the spectra of our 
PAH sample with these templates rules out Mrk231, 
NGC 1068 and Arp 220 as good matches, but the spec- 
tra are similar to those of M82 and NGC 6240. Similar 
results were found by L05, who detected strong PAH fea- 
tures in the spectra of two luminous SMGs at z ~ 2.8 
that were well fitted by an M82-type spectrum. 

3.1. The Composite SMG spectrum 

We take advantage of the similarity between the spec- 
tra in the PAH sample and of our precisely known red- 
shifts (C05) to double our S/N constructing a composite 
spectrum by averaging the individual spectra (Fig. 1). 
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Fig. 2. — Relative strengths of the 7.7-fim PAH feature as 
measured by the PAH-to-IR (8 — 1000 /im) luminosity ratio and 
the restframe EW for the SMGs in the PAH sample (large tri- 
angle s). We estimate the values for M82 {Forstcr Sc hreiber et al.l 
12003 ). Arp220 (SO), NGC6240 (A06a) and for the L05 pair of 
SMGs (small triangles) directly from their published spectra as- 
suming a linear continuum slope; the error bars reflect the uncer- 
tainties of this approach. We derive a lower EW for NGC 6240 
than A06a, due to differences in continuum definitio n. Error bars 
for th e Y05 sources, for the SB-dominated UGC5101 (Armus ct al. 
120041) and for the avera ge of 13 nearby SB galaxies studied with 
IRS {Brandl °t al - 2006) reflect the uncertainties presented by the 
authors of the respective papers. 

We use the composite spectrum to make a preliminary 
assessment of the independent contributions of SB and 
AGN activity in our PAH sample. Normalizing the lo- 
cal templates to the 7.7-fim peak in the composite spec- 
trum, we find that the composite spectrum is well fitted 
at A < 9/im by the NGC 6240 and M82 spectra (Fig. 1). 
However, neither template provides a good fit to the com- 
posite spectrum at A > 9/im: the continuum emission 
of NGC 6240 exceeds that of the composite spectrum at 
A > 12/im, while the spectrum of M82 falls below it. No 
physically reasonable additional AGN component can be 
added to the NGC 6240 spectrum to produce a good fit 
to the composite spectrum at longer wavelengths. On 
the other hand, an M82-type spectrum plus a power-law 
continuum provides a good fit to the composite SMG 
data at all wavelengths. 

3.1.1. Starburst Component 

The 7.7-/im PAH feature is generally the most promi- 
nent in the mid-IR spectra of SB galaxies. Its strength 
relative to the continuum, measured by the equivalent 
width (EW), can be used to evaluate the fractional SB 
contribution to the total bolometric output, as the hot 
mid-IR continuum is enhanced significantly in the pres- 
ence of an AGN. SB-dominated objects, such as M82 and 
NGC 6240, are characterized by larger PAH EWs than 
objects with a prominent AGN, such as Mrk231. 

EWs are sensitive to how the continuum is defined. We 
define a linear continuum by interpolating between two 
points clear of PAH emission, at 6.8 /^m and 13.7/itm, or 
at 9 fim when the spectrum does not include one of these 
points. In Fig. 2 we plot the 7.7-fim restframe EWs and 
PAH-to-IR luminosity ratios for the SMGs in the PAH 
sample with 7.7-fim coverage and for the composite spec- 
trum. The error in I^^m/Lrj^ for our SMG sample is 
dominated by a ~ 20% uncertainty in the IR luminosi- 
ties (C05) 11 . We compare a number of low- and high-z 



11 IKovacs eTaT]ll2006T) show that SMGs fall below the local FIR- 
radio relation and thus the C05-derived ijp^ values, which rely on 
this relation, are in average overestimated by a factor of ~ 2. 



sources, including 2 ULIRGs at z ~ 2 with clear PAH de- 
tections from the Y05 sample with S24^m ^ 0.9 mJy, and 
2 SMGs (SMMJ02399-0136 with S 850AI m =23mJy and 
MMJ154127+6616 with S 850AI m =14.6mJy) at z ~ 2.8 
(L05). 

According to t he line-to-con t inuum (1/c) diagnos- 
tic presented by iGenzel et al.l (|1998l ). systems with 
(l/c)7.7 Atm > 1 are classified as SB— dominated and 
those with (1/c) 7 7^ < 1, as AGN-dominated. With 
(l/ c )7-7Mni ^ 1) SBs appear to dominate the Y05, L05 and 
our PAH sample 12 . However, the distribution in 7.7-/xm 
EW and PAH-to-IR luminosity ratio in Fig. 2 may sug- 
gest a distinction in the relative SB-to-AGN contribu- 
tions, with lower values of these parameters indicating a 
stronger AGN contribution. Wc distinguish three regions 
in Fig. 2: (1) a region with l ow PAH-to-IR lumin osity ra- 
tios, occupied by Mrkl014 (| Armus et al.ll2004D and the 
24/im-bright sample of Y05; (2) an intermediate PAH- 
to-IR luminosity region where NGC 6240 and the bulk of 
the SMGs in our sample are located; and (3) a region 
with the highest PAH-to-IR luminosity ratios, occupied 
by M82 and the two SMGs in L05. 

At z ~ 2, 24-/zm flux traces 8-/zm restframe contin- 
uum; a stronger hot mid-IR continuum (produced by 
an AGN) dilutes the strength of PAH features, leading 
to lower £7.7^111/ Ljji- The location of the Y05 sample 
in the plot could follow from the selection of 24 /im- 
bright targets (S24 M , n > 0.9 mJy) at this redshift, which 
would select objects with lower SB-to-AGN ratios. SMGs 
in our sample have higher L^i^m/ ratios, similar 
to NGC 6240, which we interpret as an indication of a 
markedly stronger SB contribution to the total luminos- 
ity than the Y05 sample. With similar IR- luminosities, 
the location of the L05 SMG pair in this plot indicates 
that the 7.7-/im PAH feature is very strong. With large 
values for both the EW and the PAH-to-IR luminosity 
ratio, MM J154127 has been suggested to be dominated 
by SB-activity (L05). The large PAH-to-IR luminosity 
ratio for SMMJ02399 suggests strong SB-activity; how- 
ever, the relatively low EW value, together with the ev- 
ident strong mid-IR continuum (see Fig. 1, L05) is con- 
sistent with this source having roughly equal AGN and 
SB contributions. 

The SMGs in our PAH sample have values of EW7.7 Aim 
and L7.7 Aim /Ljr^ that place their SB-to-AGN ratio be- 
tween that of the AGN-dominated ULIRG Mrk 1014 and 
the SB M82. This is qualitatively similar to NGC 6240, 
which has both SB and AGN components. As a caveat, 
we note that even though the 7.7-/zm PAH-to-IR lumi- 
nosity ratio is associated to the SB-to-AGN ratio, it is 
also sensitive to details of the spectral energy distribu- 
tion of the system, such as the presence of multiple dust 
components at different temperatures and the amount 
of extinction (e.g. Arp220; S04). This may explain the 
particularly high 7.7-/im PAH-to-IR luminosity ratio for 
SMMJ02399. 

3.1.2. The AGN Component 

Mid-IR line diagnostics suggest that SMGs are SB- 
like. The spectra of the SB-dominated galaxies M82 and 
NGC 6240 provide a good fit to the composite spectrum 
at A < 9/im; however, only an M82-type spectrum with 

12 For our sample, (1/c) ~ 1 corresponds to EQW ~ 0.5 /an 



4 



MENENDEZ-DELMESTRE ET AL. 2006 



an additional power law AGN component gives a good 
fit to the composite spectrum at all wavelengths. The 
power-law component is defined as S v ~ v~~ 2 9 , consistent 
with the range of IR spect ral indices for 3C quasars in 
ISimpson fc Rawlingsl (|2000D . 

From the power-law component flux at 10.5 /im we esti- 
mate the X-ray luminosity (Lx ) using the correlation be- 
tween 5io.5/jm and S^-iokeV presented bv lKrabbe et al.l 
(f200lh . This yields L x ~ 10 44 erg s~ x for an AGN at the 
average redshift for the SMGs presented in this paper, 
z ~ 1.4, in reasonable agreement with the X-ray lumi- 
nosities found for the SMGs in A05. Following A05's 
approach, we compare the average X-ray-to-far-IR ratio 
of the SMGs in the PAH sample with the typical ratio for 
quasars and find that the residual flux is consistent with 
an underlying AGN contributing in the order ~10% to 
the total far-IR emission. This agrees with the A05 re- 
sult that AGN activity is often present in SMGs but does 
not dominate the energetics. Since the 10.5-/im excess is 
dominated by lower redshift sources, further SMGs at 
z < 2, included in Menendez-Delmestre ct al. (in prep), 
will better constrain this excess. 

4. CONCLUSIONS 

We present first results of a Spitzer program to char- 
acterize the mid-IR spectra of high redshift SMGs. We 
compare the spectra to well-studied local templates and 
find that SMGs have starburst mid-IR spectra more like 
M82 than the often quoted local analog Arp220. The 
composite spectrum of the SMGs in the PAH sample is 
well fitted by an M82-like starburst-component with a 
power-law continuum most likely representing a fainter 
underlying AGN. This similarity to the M82 spectrum 
suggests that the chemistry of the interstellar medium 
and radiation fields in these systems may be understood 
by looking at local galaxies in detail. Analysis of the 
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7.7- fim equivalent widths and PAH-to-IR luminosity ra- 
tios show that SMGs are markedly different from 24 /im- 
selected samples, such as the z ~ 2 sample of Y05, which 
show stronger AGN contributions. This work provides 
further evidence that SMGs host both star-formation 
and AGN activity, but that star-formation dominates the 
bolometric luminosity, reiterating the role of SMGs as 
the build-up sites for a significant fraction of the stellar 
content we see today. 

By probing the lower-redshift end of the C05 SMG 
sample distribution, this sample provides rest-frame 
wavelength coverage longwards of 9 /im to assess the 
AGN contribution. The SMG at z = 2.38, with a 
redshift closer to that of a typical C05 SMG, displays 
a potentially more AGN-dominated Mrk231-type broad 
feature at restframe ~8 /im. The difference in AGN con- 
tributions within our preliminary sample suggests an in- 
creasing relative AGN-activity in SMGs at higher red- 
shifts, potentially due to the 24/im flux limit applied to 
the sample. The full sample, with a more extended red- 
shift distribution, will provide us with additional valu- 
able information concerning the typical SMG population. 
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